Scattering of charge carriers in graphene induced by topological defects 
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We study the scattering of graphene quasiparticles by topological defects, represented by holes, 
pentagons and heptagons. For the case of holes, we obtain the phase shift and found that at low 
concentration they appear to be irrelevant for the electron transport, giving a negligible contribution 
to the resistivity. Whenever pentagons are introduced into the lattice and the fermionic current is 
constrained to move near one of them we realize that such a current is scattered with an angle that 
depends only on the number of pentagons and on the side the current taken. Such a deviation may 
be determined by means of a Young-type experiment, through the interference pattern between the 
two current branches scattered by a pentagon. In the case of a heptagon such a current is also 
scattered but it diverges from the defect, preventing a interference between two beams of current 
for the same heptagon. 

PACS numbers: 81.05.Uw, 73.63.Fg, 04.20.-q 



I. INTRODUCTION AND MOTIVATION 

Graphene is a flat monolayer of carbon atoms tightly 
packed into a two-dimensional {2D) honeycomb lattice, 
(Fig.[T|) or it can be viewed as an individual atomic plane 
pulled out of bulk graphiteii^. It is the first example of 
a truly atomic two-dimensional {2D) crystalline system 
and it is the basic building blocks for graphitic materi- 
als such as fullerenes (graphene balled into a sphere) or 
carbon nanotubes (graphene roUed-up in cylinders )^^—. 
Additionally, it is a zero-gap semiconductor, in which 
the low energy spectrum is correctly described by the 
(2 + 1)D Dirac-like equation for a massless particle^i^i^ 
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where is the Fermi velocity, which plays the role of 
the speed of light {vp « c/300), a = {ax,o-y) are the 
2D Pauli matrices, p ~ —ih\7 is the linear momentum 
operator and |^') is a two-component spinor. Therefore, 
the quasiparticles can be viewed as electrons that have 
lost their masses or as (massless) neutrinos that acquired 
the electronic charge. Such a spectrum makes graphene 
a material with unique electronic properties. Its descrip- 
tion by means of the Dirac equation is a direct conse- 
quence of graphene?s crystal symmetry. Its honeycomb 
lattice is made up of two equivalent triangular carbon 
sublattices A and B, whereas its cosine- like energy bands 
associated with the sublattices intersect at zero energy 
{E = 0) near the edges of the Brillouin zone, giving rise to 
a conical section spectrum at low energies^, say, |£'| < 1 
eV. In this honeycomb lattice, the two-component spinor 

4') is referred to as pseudospin since it is an index in- 
dicating two interpenetrating triangular sublattices A e 
B, which is similar to spin index (up and down) in quan- 
tum electrodynamics (QED). It is common to regard the 
sublattice degree of freedom as a pseudospin, with the 
A sublattice being the "up", |-|-) and B sublattice being 
the "down", |— ), i.e.. 
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FIG. 1: (Color online) Graphene is a flat monolayer of carbon 
atoms having a honeycomb lattice, consisting of two interpen- 
etrating triangular sublattices (green and red circles) . Exper- 
imental techniques provide high-quality graphene crystallites 
up to 100/im in size, which is sufficient for most research pur- 
pose, including the ones considered here. In this "perfect" 
layer the charge carries can travel thousands of interatomic 
distances without scattering. 
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Since vp <^ c, it is a slow relativistic system or a strong 
coupling version of QED since the graphene's dimen- 
sionless coupling constant, e^/hvp w 1 much higher 
than its QED analogue, the fine structure constant 
e^/hc w 1/137. All these properties make the graphene 
a very interesting system, which provides a way to probe 
QED phenomena, for instance, by measuring its elec- 
tronic properties. Several proposals for testing some 
predicted, but not yet observed phenomena in QED, in- 
cluding the Klein parados^, vacuum polarizatiork^ and 
atomic coUapscji^, arc some topics under investigation in 
graphene. 

Here, we would like to study the behavior of graphene 
quasiparticles in the presence of defects in the crystalline 
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structure of the material. We shall consider three types 
of defects: holes, pentagons and heptagons. All these de- 
fects can be incorporated by removing or inserting a few 
carbon atoms. Understanding how these defects modify 
the transport properties of graphene is crucial to achieve 
future electronic devices using carbon-made materials. 

The presence of defects like pentagons (heptagons) in- 
duces positive (negative) curvature in the material. At 
some extent the charge carriers motion in the presence of 
pentagonal (heptagonal) defects is identical to fermions 
moving in a (2 + 1)D gravitational space-time generated 
by positive ("negative") point-like masses. Then we may 
employ results from gravity to analyze some effects con- 
cerning charge carriers in the presence of such defects. 
For instance, in the presence of pentagons and heptagons 
we found that the phase shifts of the carriers moving near 
these defects depend only on the number of pentagons 
or heptagons in the graphene sheet and are identical for 
charge carriers in the sublattice A or B. In the case of 
the holes we computed the phase shifts for the scattered 
electrons and found that at low concentration they are ir- 
relevant for the electrons transport giving no appreciable 
contribution to the resistivity. 



3-dimensional Levi-Civita symbol (e"-'^^ = -|-1). (The 
word covariant must be used carefully because vp is 
not invariant, being only a parameter and the term 
covariant is used to refer only the form the equation 
is written.) We may expand the solutions of the free 
massless Dirac equation in plane waves, once rotational 
invariancc allows to separate the 9 variable, so that the 
diagonalized angular momentum = —ih-^ + f'^'^j 
yielding partial waves with angular momentum {n+ ^)h, 
takes the form, 
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The components of the radial spinor Un(r), given by /«(?") 
and gn{f)i satisfy the Besscl equation of order n and n+1, 
say: 



dr^ r dr 
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II. CHARGE CARRIERS DYNAMICS IN A 
NON-SIMPLY CONNECTED GRAPHENE SHEET 

In this section we discuss the scattering of the charge 
carriers in a non-simply connected graphene sheet. In 
the continuum model for the graphene, it is assumed that 
there is a hole of radius rg cut from the system center, lo- 
cated at the origin. Thus the motion of the quasiparticles 
is performed on a flat 2D support given by a non-simply 
connected manifold, which can be viewed as defect in the 
material. This model allows for the investigation of scat- 
tering effects as a function of the hole radius and could 
shed some light on the high charge carrier mobilities, fact 
observed in graphene^iiii. In addition, it is interesting to 
compare our results with the scattering by a short-range 
potential because it appears that in the case of graphene 
the contribution of such small defects to the resistivity is 
essentially smaller than for conventional nonrelativistic 
two-dimensional electron gas^iii. 

To determine the phase-shift of the scattered wave 
function as well as the scattering cross section one has 
to solve the two-dimensional Dirac equation ([T|) which, 
for the case of massless particles, can be writing in a 
covariant form, like below: 

ih-i^d^iij{x) = 0, (3) 



where n = , ±1 , ±2 , . . . is the angular-momentum 
number and k — > (we arc considering only solu- 
tions with E > Q that describe the electronic dynamics). 
Remember that the spinor index in graphene labels its 
two crystal sublattices rather than directions of the ac- 
tual spin. The solution for the radial spinor outside of 
the hole (i.e., for r > rg) is given by: 



( Bin Jn{kr) + B2n Njkr) \ 
\ Bsn Jn+i (kr) + Bin iV„+i (fcr) ^ ' 

where J„ and Nn are the Besscl functions of first and 
second kinds (Neumann function), respectively, whereas 
Bjn {j = 1 J 2 , 3 , 4) are constants. These constants 
must be determined by the appropriate boundary con- 
ditions specified to completely define the problem. From 
the physical point of view, the correct boundary condi- 
tion is determined by the requirement of vanishing the 
net energy flux into the hole, which is a region absent of 
lattice degrees of freedom. Consequently, the fields must 
arrange themselves in such a way that the energy flux 
from the incoming modes (asymptotically behaving like 
g-ifcr-j exactly cancel that from the outgoing waves. 
Let us recall that in the case without a hole, one has: 



where the covariant derivative is = 
[{l/vF)d/dt , d/dx , d/dy], the 7-matrices are 
7*^ ~ cr'^, 7"'^ = i(7^ and 7'^ = —ia^, obeying 
7^'^'^ = jyM!^ _ je^''"7a, if is the Minkowski ten- 
sor metric, diag(?7^'^) — (+1, —1, —1) and e'^"" is the 



To determine the phase-shift (5„, we need the asymptotic 
behavior, r 00, of the radial spinor ^ sa.y^: 
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where the spinor indices A and B labels the two crys- 
tal sublattices. Now, imposing Neumann boundary con- 



dition (NBC) on the wave-functions, 
obtain: 



0, we 



Fermi level, taking into account the four-fold degeneracy, 
the spin and valley of the graphene, and r is the mean- 
free path. For electrons of the sublattice A we have: 



da{0) 
d9 
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The Dirac eq. ([3]) or ([5]) has an important symmetry 
under the interchanging / i — > g , n < — > —n — 1 which 
implies that i„ = t-n-i- Thus, eq. ([TO)) can be written 
in the form: 



Unir) = 



where 



Si„[J„(fcr) - tan(in(fcr-o))A^„(fcr)] 
-B3„[j„+i(fcr) - tan{tn+iikro))Nn+i{kr)] 



tan[t„(fcro)] 



(10) 



(11) 



The terms proportional to Bessel (Neumann) functions 
describe the incident (scattered) waves. Comparing their 
asymptotic behavior with those for plane waves, eq. ([9]), 
we obtain the phase shifts. Asymptotically, the scattered 
wave- function M„(r) behaves like follows: 



d<j{6) _ 2 

de ~ ^ 
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(17) 



For small energies, kro 1, which is typical for 
graphene, one has 



[(2ri + l)!!]2(2n + 3)(n + l) 



(18) 



and thus the s-scattering {n = 0) dominates. With eqs. 
(jlTHlSp the impurities contribution to the resistivity may 
be estimated as 



2 ( n-K TT , , 

— — COS fcr + <„(fcro, 

Tr/cr \ 2 4 
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+i„+i(fcro)) , (12) 



giving us the phase-shift 5n of the n-th partial wave that 
completely determines the fermionic scattering: 



tn{kro) , 



tn+i{kro) . 



(13) 
(14) 



Now, if we consider a small concentration of point-like 
defects with concentration n^cf angie-depedent scattering 
cross section, (t(0), their contribution to resistivity, p, 
reads^: 



P = 



T{k[ 



e^vlN{EF)T{kF) ' 

ndefVF de^^{l^cos9), (15) 







d9 



where N{Ef) — 2kF/iThvF is the density of states at the 



P - ndei 



e^kki 



-{krof 



(19) 



This means that the scattering induced by small holes 
(with radius around a few angstroms; some lattice spac- 
ings) at low concentration are irrelevant for the electronic 
transport in graphene, giving a negligible contribution to 
the resistivity. For the case of a potential V{r) = Vb at 
r < Rq and V{r) = a.t r > Rq, the estimation for 
this type of impurity contribution to the resistivity is^ 
p ~ (/i/4e^)ndcf^0' gi'^ing ^ negligible contribution to 
the resistivity when the radius of the potential Rq is of 
the order of interatomic distances and at a low concen- 
tration as above. 

For an intuitive understanding of the result (fT9|l let 
us recall that light does not experience obstacles with 
sizes much smaller than its wavelength, whereas massless 
Dirac electrons have the same dispersion relation behav- 
ing as light in some aspects (duality wave-particle). This 
same interpretation can explain the results for a short 
range potentialii. Ours results for the resistivity are in 
agreement with those presented in Ref.[ll|, where the 
authors show that intrinsic corrugations of a graphene 
sheet create a long-range scattering potential and lead to 
significant resistivity that could explain the existing ex- 
perimental data about the graphene resistivity (to more 
details see Ref.[ll|). For more details about scattering 
of the charge carriers by defects in graphene see, for in- 
stance [Hil-lill and references therein. 
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III. CHARGE-CARRIERS DYNAMICS IN THE 
PRESENCE OF PENTAGONAL AND 
HEPTAGONAL DEFECTS 



ever building the cones, say, the number of inserted pen- 
tagons must be kept a minimum. 



Here, pentagons and heptagons defects are considered 
and their effects in the charge carriers dynamics is in- 
vestigated. Substitution of an hexagon by other type of 
polygon with n = 6 — Ud sides, where is an inte- 
ger smaller than 6, in the lattice without affecting the 
threefold coordination of the carbon atoms leades to the 
wrapping of the graphene sheet. These defects can be 
seen as disclinations of the lattice locally acquiring a fi- 
nite curvature. The accumulations of various defects may 
lead to closed shapes as fullerenes. Rings with n < 6 sides 
(ud > 0) give rise to positively curved structures whereas 
polygons with n > 6 sides {ud < 0) yields negative curva- 
ture. This induced curvature holds only near the defect, 
whereas the graphene sheet is flat away from the defect 
itself, as the conical surface is out from the ape3*i^"— . 

Exploring the 2D character and flexibility of this ma- 
terial, our idea is to propose a system which one or more 
sectors are excised from a graphene and the remainder is 
joined seamlessly, (Fig. [2]). In fact, the missed link of each 
carbon atom resting at the two edges of the remaining 
graphene sheet can be, in principle, covalently bounded. 
The nucleation and growth of curved carbon structures 
remain to be well-understood. It is claimed that the oc- 
currence of pentagons, which yields 60° disclination de- 
fects in hexagonal graphitic network, is a key element 
in this scenario. Particularly, considering the symmetry 
of a graphitic sheet and the Euler theorem, it can be 
shown that only five types of cones (incorporating one 
to five pentagons) can be made of a continuous graphene 
shee t^^'^° . In the case of a cone with Ud > 0, the value 
nd, [nd = 1 , . . . , 5) is related to the conical angle 7: 



■ 1 . rid 
sm — = i — . 



(20) 



The deficit-angle induced by the conical singularity is 
given by 27r(l — nd/6)- The pentagonal defect can be 
presented as a pseudo-magnetic vortex at the apex of a 
graphitic cone, being the flux of the vortex related to 
the deficit angle of the cone (see RefJ^). The possi- 
ble five graphitic cones mentioned earlier are then given 
by 7 = 19.2°, 38.9°, 60°, 84.6°, 112.9°i^i^. Cones with a 
heptagon have negative curvature and are obtained by a 
insertion of a angular sector in the carbon sheet. Then if 
Ud < 0, —Hd counts the number of such sectors inserted 
into the graphene sheet. Our aim is, therefore, to see the 
influences that such a special graphene structure could 
induce on quasiparticlc wavefunctions (spinors); surely, 
these influences may create new perspectives in the elec- 
tronic transport properties, which are determined by the 
quasiparticles constrained to move on the conical surface. 
Once only the honeycomb lattice ensures the graphene 
peculiar dispersion relation E = vpp, lattice distortions 
deviations should be minimized as much as possible wher- 
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FIG. 2: (Color online) Removing a wedge from the graphene 
(left) and identifying the edges (right), a cone results. The 
motion of the charge carriers in an ideal conical graphene 
is equivalent to that of a massless Dirac particle in a grav- 
itational field of a static particle of mass M in a (2 -I- 1)D 
space-time. 



To study the scattering of the carriers in graphene 
by topological defects we employ the analogy between 
defects in condensed matter physics and in (2 -I- 1)- 
dimensional gravity^! as far as possible. For example, the 
dynamics of charge carriers in an ideal conical graphene is 
equivalent to that of a massless Dirac particle in a grav- 
itational field of a static point-like mass in a (2 -f- 1)D 
space-timei^i^i^. Specifically, we shall consider one of 
the simplest curved manifold, which is associated to the 
Schwarzschild solution in (2+1) dimensions: a space-time 
locally flat with global nontrivial properties described be- 
low. 

To proceed further, we summarize some aspects of gen- 
eral relativity in three space-time dimensions. It funda- 
mentally differs from its 4-dimensional counterpart. In- 
deed, it exhibits some unusual features, which can be 
deduced from the properties of the Einstein field equa- 
tions and the Riemann curvature tensor i?^^^^— In 
regions free of matter (where the momentum-energy ten- 
sor T**^ vanishes), the space-time is locally flat when the 
cosmological constant is zero (the Einstein tensor G^'' 
vanishes). However, this does not mean that a massive 
source has no gravitational effects: a light beam passing 
by a massive, point-like mass will be deflected^ii"— and 
parallel transport in a closed circuit around it will in gen- 
erally gives nontrivial results^^i^. Indeed, while the local 
curvature vanishes outside the sources, there are nontriv- 
ial global effects. For instance, for the special case of a 
point-like mass, M, sitting at rest in the origin, the line 
element is given by 



ds^ = - dp^ - p^a^d9^ 



(21) 



with < 6* < 27r and a = 1 - 4GM (G is the New- 
ton's constant, with dimensions of [length]^ instead of 



[length] , in natural units h 



1). Note that al- 



though the situation looks trivial, the coordinate 9 ranges 
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FIG. 3: (Color online) Classical trajectories (bue and red) of 
the scattered particles deflecteds by heptagons (above) and 
pentagons (below). The trajectories are the asymptotic mo- 
tion on cone projected onto x — y plane of the embedding 
three dimensional space, see the text. 



a digression concerning scattering of the charge carriers in 
graphene as classical relativistic particles. The classical 
equation of motion, determined by relativistic geodesic 
equation for the particles in a cone reads i + F'' 



0, where the overdot indicates differentiation with respect 
to any convenient affine variable r that parametrizes the 
path a;^(r)^. The angle of scattering ±a; for the motion 
of the particles in a cone can be obtained by integration 
of the classical equations of motion and is given by2^. 



±w = ±7r(a-^ - I) = ±7r 



/3 



1-/3' 



(23) 



where ± refers to the side the charge carriers (current) 
trajectory pass around the defect, see (Fig. Note 
that the result above is valid for all values of /3 despite 
the sector was removed or inserted. The scattering an- 
gle above, presented in the embedded coordinate system 
measures the deflection of the asymptotic motion on 
cone projected onto x — y plane of the embedding three 
dimensional space. The result above suggests that a pen- 
tagon or heptagon may be used for deviating the planar 
current in graphene. 

To obtain the correct current deviations in graphene we 
have to solve the Dirac equation ([T]) defined in a cone, 



from to 27rQ!, indicating an angular deficit in space. 
Then, the spatial part of the metric is that of a plane with 
a wedge removed and edges identified; the unique 2D spa- 
tial geometry satisfying this description is the cone^i. 

Alternatively, one may use embedded coordinates r 
and 9 in the three-dimensional Euclidian space which 
extend over the complete range, < r < oo , < 
< 27r, and describe a cone with the constraint z = 
^ (a"2 — l)(a:;2 + 1/^), being the line element given by2^: 

ds^ = di^ - oT"^ dr^ - dd^ . (22) 

The attributes of the source are coded in the global prop- 
erties of the locally flat variables. All the information 
lies in the non-trivial boundary conditions, which is im- 
portant for the quantum scattering of graphene charge 
carriers by defects, like pentagons and heptagons. 

In the case of a topological defect in graphene, it is use- 
ful to change the gravitational term AGM by the symbol 
/3, so that 2tt/3 (for < /3 < 1) gives the deficit of angle 
measuring the magnitude of the removed sector whereas 
— 27r/3 (for — oo < f3 < 0) accounts for the angle in excess 
associated to the insertion of a sector. The parameter f3 
takes only discrete values because of the lattice symmetry 
of the graphene as discussed after cq. (|20|) . 

Before analyzing the quantum mechanical scattering 
by conical defects (pentagons and heptagons) let us make 



ih^'^Ea "D^^jj = , (24) 

where D,, = 9^ + icJ^;afc^T"^ d'^'' = i[7°,7'], and Ea" 
is the dreibein in coordinates {t,r,9). The spin connec- 
tion ij^;ab = —'^^l■.ba may be written in three dimensions 
as cijpiab — iabci^n with eabc the Levi-Civita symbol 
as bcforo2^. The rotational invariance of the problem 
enables us to choose positive energy solutions that are 
simultaneously angular momentum eigenfunctions, with 
eigenvalue {n + ^)h: 

u„{r)e-^^'/'^ = e'(»+i-i-')'^ ( ^ ) e"'^*/'' , (25) 
where n — 0, ±1, ±2, .... The solutions for E' > are^^: 
u^{r)^{e^rM^r), (26) 
<{r) = (e„)"+V.(«0- (27) 

Here, is the Bessel function of order v = e„/a(n+ (1 =F 
a)/2), n = 0, ±1, ±2,. . ., K = E/hvpa.E > 0,e„ = ±1 
and the same sign has to be chosen for the upper and 
lower components of the spinor Un{r). For < a < 1 or 
< /3 < 1 (remember that a = 1 — /?) we must choose 
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e„ = sign(n+(l— a)/2) = sign n (sign = 1) to have both 
components regular at the origin. The asymptotic form 
of the Bessel functions determines the phase shifts (they 
are identical for the upper and lower components) and are 
given by2^: 



the pseudospin, may suffer important changes in the con- 
ical graphene because a hopping of a quasiparticle, which 
was previously, say, in the A sublattice, will make it to 
become out of phase with all quasiparticles occupying the 
B sublattice. 



5n = -e„ — ((l-a)n+(l-a)/2) 



e„ = sign(n + (1 - a)/2) = sign(n + /3/2) . (29) 

The phase-shifts depend only on the number of sectors 
removed or inserted in the graphene sheet, accounted by 
a = 1 — p. If — oo < /3 < 0, we need to be careful 
because e„ = ±1 depending upon the value of (n -I- /3/2) 
(but the phase-shifts remains as above) and the phase- 
shifts depends only on the number of sectors (heptagons) 
inserted in the flat graphene sheet. In the presence of 
heptagons the carriers dynamics is identical to the that 
movement of the electrons in the gravitational field of a 
negative mass (although not possible in gravitation, this 
is feasible in the present context). 

Note that the phase shift ([28]) measures the deflec- 
tion of the asymptotic motion on the cone projected onto 
X — y plane being qualitatively identical to the classical 
case discussed before. When there is a pentagon into the 
lattice and the fermionic current is constrained to pass 
around and sufficiently close to it such a current is scat- 
tered by the defect with an angle which depends only on 
the number of sectors removed in the graphene and on 
the side current passed (see Fig. [3]). After passing by 
the pentagon the scattered current trajectories cross and 
yields an interference pattern. In the case of a heptagon, 
such a current is scattered but the trajectories diverge 
each other. 

The presence of pentagons or heptagons in the sheet of 
graphene may manifest as fluctuations in the concentra- 
tion of charge carriers, modifying several physical proper- 
ties. For instance, in a planar graphene, it is known that 
the hopping of electrons between sublattices produces an 
effective magnetic field which is proportional, in magni- 
tude and direction, to the momentum measured from the 
Brillouin-zone corners. This effective field, which acts on 



IV. CONCLUSIONS 

We have studied the scattering of graphene quasiparti- 
cles by topological defects like holes, pentagons and hep- 
tagons. Wc obtain the phase shift of the wave-function 
in all cases. For the case of holes, the main contribu- 
tion concerns the s scattering and even in this case they 
do not change the resistivity of the sample, at least at 
low concentrations (like occurs to short range potential 
impurities) . 

Wc realize that when the fermionic current is con- 
strained to move near and around of pentagons and hep- 
tagons introduced in the lattice, it is scattered with an 
angle that depends only on the number of defects and 
on which side the current taken. Such a deviation may 
be determined by means of a Young- type experiment, 
through the interference pattern between the two cur- 
rents scattered by the pentagon. In the case of a hep- 
tagon such a current is also scattered but it diverges from 
the defect. Perhaps these effects could be used to build 
channels of currents in future electronic devices using 
graphene. In addition, graphene would provide an ap- 
pealing way to experimentally explore general relativity 
in two spatial dimensions since such effects are predicted 
by this tlieoryi^>^>2a. 

Conical spaces generated by topological defects like 
vortices have been recently considered and the quan- 
tum mechanical scattering of a nonrelativistic particle- 
have been considered in the work of Ref.— with results 
qualitatively similar to those obtained here. 
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